Self-trapping threshold in disordered nonlinear photonic lattices 
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We investigate numerically and experimentally the influence of coupling disorder on the self-trapping dynamics in nonlinear 
one-dimensional optical waveguide arrays. The existence of a lower and upper bound of the effective average propagation 
constant allows for a generalized definition of the threshold power for the onset of soliton localization. When compared to 
perfectly ordered systems, this threshold is found to decrease in the presence of coupling disorder. © 2013 Optical Society 
of America 
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Discrete solitons \\^ emerge in periodic nonlinear lattices 
when self-focusing becomes sufficiently strong to balance 
diffraction. In spite of this very general mechanism, the transi- 
tion dynamics between extended and localized states depend 
strongly on the geometry and dimensionality of the specific 
system. For example, the critical self-trapping nonlinearity 
in a chain of coupled oscillators decreases with the length 
of the chain, as can seen comparing the values of a coupled 
dimerpl, smaller one-dimensional chains ||4| and very long 
ones 13]. A first attempt to find a more general criterion led to 
the discovery of a common ratio of^ 1.3 between the critical 
nonlinearity and the minimum bound- state energy for defect 
states in lattices with different dimensions and topologies 1 6l . 
Experimentally, the transition was observed in arrays of Al- 
GaAs 1 7 1, in arrays made with the fs-laser direct- writing tech- 
nology |8| and nematic liquid crystal arrays f9l. 
Disorder can manifest itself as random change to the lo- 
cal properties of all lattice sites. Known as Anderson local- 
ization, the ensuing interference between multiple scattering 
paths may significantly impede and even suppress the trans- 
verse transport of waves |T0| . Exponential localization was 
first directly observed in disordered photonic lattices fll||12|. 
Similarly, random modifications of the coupling between lat- 
tice sites can lead to exponentially localized states |T3| , as 
was demonstrated in waveguide arrays 1 14] where coupling 
disorder corresponds to variations of the waveguide separa- 
tion (TSj. It is known that Anderson localization can be de- 
stroyed by weak nonlinearity, giving rise to unlimited subdif- 
fusive spreading for very large systems [16, - 1 8 J . On the other 
hand, when the system size is smaller than the (linear) local- 
ization length, the interplay of diagonal disorder and nonlin- 
earity induces a smoothening in the distribution of a spread- 
ing initial single- site excitation p9| . In the limit of large non- 
linearities, self- trapping of a large fraction of the initial wave 
packet occurs irrespective of disorder [20,21 1. 
The question arises as to how the simultaneous action of both 
nonlinearity and disorder influences not only the width of 
the output |[TT][T2][22|, but especially the dynamics of self- 
trapping. The existence of a lower and upper bound of the 



average propagation constant was shown recently for various 
ordered lattice configurations of different dimensions, leading 
to the definition of a generalized power threshold |23|. In this 
letter, we show numerically and verify experimentally that the 
presence of coupling disorder entails a systematic decrease of 
this threshold for dynamical excitation in one-dimensional op- 
tical waveguide arrays. 

In the coupled-mode approximation, the evolution of light in a 
ID lattice of N waveguides along the propagation direction z 
is described by a discrete nonlinear Schrodinger-like equation: 
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where Un corresponds to the modal field amplitude at site n, 
Cn is the coupling between guides n and ^ + 1, and 7=1 
(7 = — 1) is the nonlinear coefficient for focusing (defocusing) 
nonlinearity. Hereafter, we take the nonlinear coefficient 7=1 
and consider primarily the case without diagonal disorder, thus 
e„ = . The effective strength of the nonlinearity is then given 
by the conserved total power (norm) P = Y.nWA'^ • Further- 
more, we will consider single-site excitations, i.e. Un{z = 0) = 
VPSfi^no , where no corresponds to a central position in the lat- 
tice. The second conserved quantity of the system, the Hamil- 
tonian 
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then has the value of H = P^/2. In order to evaluate the size 
of the wave packet, we use the participation number 



R: 
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which approaches 1 for a highly localized wave packet and 
tends to N for an extended excitation. R serves as a measure 
of how many lattice sites effectively contribute to a given pro- 
file. In comparison to other quantities such as the variance or 
standard deviation, the participation number emphasizes sub- 
stantially excited sites and is less sensitive to weak excitations 
that may propagate far away from the initial site. 




Fig. 1. Power-dependent effective propagation constant. 
Shown is the mean (A^) (solid black graph) and A^ \z=n/a (solid 
gray) obtained from central single- site excitations in a homo- 
geneous lattice (Cn = 1) with N = 221. Inset: {R) vs. P. 

In analogy to stationary solutions Un{z) = w^exp (/Az), we as- 
sign an effective propagation constant A^ to any momentary 
profile. This quantity can be obtained by multiplying Eq.([T]) 
by w* and taking the sum over all lattice sites: 

n n 

This expression is related to the Hamiltonian: 

2^' ' 2 R 

After inserting H = P^/2, we finally arrive at 



2^ + ^ 



(4) 



Let us first consider the limiting cases. In the strongly non- 
linear regime, diffraction is totally suppressed and the wave 
packet remains localized with a participation ratio of R ^ 
1; the effective propagation constant becomes equal to the 
power: A^ = P. Conversely, in a diffractive process, the wave 
packet spreads across the whole lattice. The value of R thus 
increases up to a value of the order of N, yielding A^ = P/2 
for sufficiently large systems. 

Figure [T] illustrates the power-dependent behavior of A^ for a 
single-site initial excitation at the center of a uniform lattice 
with N = 22\ waveguides, the inset shows {R) vs. P. At van- 
ishing excitation power, the 5 -like input pattern corresponds 
to a flat excitation of all modes of the lattice; A^ thus consti- 
tutes the average over the propagation band and emerges at 
its center (A^ = 0). With increasing power, nonlinear contri- 
butions become relevant as the modes start to interact and to 
exchange energy. Consequently, A^ exhibits oscillations along 
z with a power-dependent period. It has been shown fST] that 
Xe represents the average propagation constant of all excited 
modes. Its behavior inside the sector P/2 ^ A^ ^ P serves as 
indicator as to whether the wave packet's evolution is dom- 
inated by diffraction or localization. The power-dependent 



transition between those two regimes becomes clearly visi- 
ble when the rapid oscillations of A^ {P) are removed by av- 
eraging over a certain interval of propagation distances (here: 
z G [N/10,N /A]). In the following, we define the threshold 
power Pth as the lowest value where A^ exceeds the linear limit 
of P/2 by a certain fixed cutoff value A. 
We now turn to lattices with coupling disorder, where the in- 
dividual coupling constants Q = 1 + /i are uniformly distrib- 
uted with ji e [-Wc/2,Wc/2]. Sets of 50 realizations each for 
different degrees of disorder were independently evaluated. 
Figure [2ja) shows the ensemble-averaged value of (A^) —P/2 
for off-disorder. For small P, all curves remain close to zero. 
Note, however, that the most strongly disordered ensemble 
(with W = 0.6) exhibits the largest values. This illustrates the 
influence of Anderson localization, which gives rise to less 
extended linear modes. Consequently the condition R:^ I is 
fulfilled to a lesser degree, for increasing disorder. The de- 
crease for on-diagonal disorder is weaker. Due to the lack of 
a threshold for the existence of nonlinear stationary modes in 
ID lattices, the onset of localization is gradual p3| . Neverthe- 
less, a pronounced transition is visible; a choice of A = 0.05 
provides a good measure for the departure from the linear 
behavior. Clearly, disorder serves to smoothen the transition 
in the ensemble averages, but also systematically lowers the 
power threshold for the onset of nonlinear localization [see 
black curve in Fig. |2jb)]. For comparison, the threshold for 
equivalent amounts of on-diagonal disorder is shown [dashed 
curve in Fig. [2jb)]. 

To verify our results experimentally, we fabricated waveg- 
uide arrays in a 100 mm long fused silica sample using the 
fs-laser direct- writing technology |^|. The coupling disorder 




Fig. 2. (Color online) (a) Transition between linear behav- 
ior and nonlinear localization for varying off-diagonal disor- 
der. Shown is (A^) - P/2 for W^ = ( solid black), W^ = 0.1 
(dashed blue) and Wc = 0.6 (dotted red), (b) Ensemble aver- 
ages and standard error of Pth obtained for A = 0.05 for off- 
(on-)diagonal disorder with full (dashed) lines. 




(a) 



0.2 0.4 0^6 0.: 

Average power P [mW] 




(b) 



0.2 0.4 0.6 0.8 

Average power P [mW] 



1.0 



Fig. 3. (Color online) (a) Experimental mean output intensity 
distributions at z = 1.28;r/2C for a = C [top], a = 0.4 C 
[center] and O" = 0.8 C [bottom], (b) Transition region: Power- 
dependent deviation from linear behavior (A^) — P/2 for G = 
C (solid black), a = 0.4 C (dashed blue) and a = 0.8 C 
(dotted red). Inset: The effective propagation constant remains 
confined to the interval P/2 ^ A^ ^ P. 

C = C ± cr was realized for three different values ((7 = 0, 
0.4 C and 0.8C) by varying the waveguide separations |15| 
and ensembles of 18 independent realizations were exam- 
ined. The sample length corresponded to a propagation dis- 
tance of 1.287r/2C with C = 0.02 mm~^ Nonlinear excita- 
tion was achieved by a Ti: sapphire laser system, delivering 
300 fs pulses with a repetition rate of IkHz at 800 nm. The 
plots in Fig. [3ja) show the ensemble- averaged output inten- 
sity distributions obtained with average input powers between 
< P < 1.0 mW for vanishing, intermediate and high disor- 
der. Despite the diffractive background from the pulsed excita- 
tion, the accelerated onset of localization is clearly visible for 
the disordered cases. Figure[3jb) substantiates this observation 
by means of the effective propagation constant obtained from 
the patterns. In agreement with our simulations, the dynamics 
are generally confined to the region P/2 ^ A^ ^ P [see inset] . 
Note that due to the limited propagation length, A^ systemat- 
ically exceeds P/2. Nevertheless, the transition to the regime 
of nonlinear localization is clearly visible in all cases, and is 
shifted towards lower P for larger degrees of disorder a. 
In conclusion, we have analyzed the nonlinear localization dy- 
namics in disordered photonic lattices. Based on the definition 
of a dynamical threshold that relies on analytical estimates of 
the upper and lower bounds of the propagation constant, we 



have shown numerically and verified experimentally that cou- 
pling disorder may strongly inhibit transport, causing fewer 
sites to be substantially excited and systematically decreases 
the threshold power required to dynamically excite nonlinear 
localized wave packets. 
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